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I. INTRODUCTION 

Entanglement of quantum states is the most non-classical 
feature of quantum systems. It shows up when the system 
consists of two (or more) subsystems and the total state cannot 
be written as a product state. This notion can be generalized 
to mixed states, and mixed state is non-separable or entan- 
gled if a corresponding density matrix cannot be expressed 
as convex combination of tensor products of density matrices 
of subsystems Ol. f^ire entangled states as superpositions of 
multiparticle states, are fragile with respect to noise resulting 
from interaction with environment. So to control the effects 
of noise, it is important to understand details of the process 
of disentanglement i.e. to analyse how entanglement can be 
destroyed by this interaction. 

In the paper, we describe the model of two two-level atoms 
interacting with thermal bath at formally infinite temperature. 
When atoms are separated by large distance, we can assume 
that the atoms are located inside two independent baths. In 
such a case time evolution of two-atomic system is given, 
in the Markovian approximation, by the ergodic dynamical 
semi-group i.e. a trace preserving semi-group of completely 
positive operators which have the maximally mixed state ^ I4 
invariant. Density matrix corresponding to the state of the 
system satisfies the master equation which right hand side is 
given by the Lindblad generator of the semi-group. Since two 
atoms are largely separated, the dipol-dipol interaction and 
the photon exchanges between atoms are negligible, so in our 
model the generator is parametrized only by the dissipation 
rate T. 

From the general properties of this kind of evolution it fol- 
lows that the entanglement as a function of time always de- 
creases to zero, and that dynamics needs only a finite time to 
disentangle any initially entangled state. Thus for all times 
t greater then some time of disentanglement fd(p), the states 
p(f ) are separable. This finite-time effect should be compared 
with asymptotic noise-induced decoherence effects (see also 
Q IH). We calculate time evolution of some class of initial 
density matrices and obtain analytic expression for its entan- 
glement. We also find formulas for in the cases of pure ini- 
tial states and some mixed initial states. In that examples, the 
time of disentanglement is the increasing function of the initial 
entanglement. Another interesting aspect of dissipative evolu- 
tion studied in the paper is connected with nonlocal properties 
of quantum states. Nonlocality of quantum theory manifests 



by violation of Bell inequalities, and in the case of two two- 
level atoms can be quantified by some parameter ranging from 
for local states to 1 for states maximally violating some Bell 
inequality. Our dynamics enables also to consider evolution 
of this parameter In particular, we show that the time after 
which all nonlocal properties of quantum state is lost is much 
shorter then the time of disentanglement. 



II. TIME EVOLUTION IN A NOISY ENVIRONMENT 

A. Two-level atom in a noisy environment 

Time evolution of a density matrix of two-level atom A in- 
side the bath with finite temperature T can be described by the 
following master equation 
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where 



ri{[cT-,pa+] + [o_p,o+]} 



o± = - (ai ± /02) 



and we identify ground state |0) and excited state |1) of the 

Moreover, 



atom with vectors ( ^ j and | q ) in ^^(4 



where 



r|=yon((oo), =yo(l +n(coo)) 



(11.2) 



yo is a spontaneous emission rate of the atom and coo is the 
frequency of the transition |0) Since 



1 



when P 0, for very high temperature {T °°) we can as- 
sume that |[4|] 
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In this case (II. 1) reduces to 
dp 



satisfy 



— = Lrp = r (a+po_ + o_po+ - p) 



(11.3) 



and Ly generates the semi-group 

T, = e'-^r 

which is ergodic i.e. {r,}(>o has a unique asymptotic state 
which is maximally mixed state ^ in C^. So the relaxation 
process described by (II. 3) which models the evolution of two- 
level atom inside the bath with very high (infinite) temper- 
ature, brings all initial states of the atom into the state with 
maximal entropy. In the other words, the semi-group gener- 
ated by (II. 3) describes open quantum system (two level atom) 
interacting with noisy environment (a bath with very high tem- 
perature). 

B. Two independent two-level atoms 

In the case of two separated two-level atoms A and B lo- 
cated inside independent baths with formally infinite temper- 
atures, the generator of the corresponding semi-group is given 
by the following generalization of (II. 1) 



dp 



dt 



= Lfp^" =r (o^p^^c/ + p^ V+ 
o«p'^V+a«p'^X-2p 



(11.4) 



where p"^^ is the state of the compound system AB, described 
by the Hilbert space ^{ab = J^a® — and the algebra of 
observables 



aAB-M2x2(C)®M2x2(C)~M4x4(C) 



and 



o4 =a±( 
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The semi-group generated by Lp is also ergodic with unique 

asymptotic state ^ in C*. To simplify the discussion of the 
evolution of two-atomic system, let us introduce the basis of 
so called collective states in the Hilbert space ^Kab = 1 5]. 
If 
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1 1) , /2 - 1 1) §5 |0) , /3 = |0) ® 1 1) , /4 = |0) ® |0) 



then this basis containing excited state, ground state and sym- 
metric and antisymmetric combination of the product states, 
is defined as follows 

k) = /l , \g) = /4, \S) = -J= (/2 +/3), \a) = -J= (/2 -/3) 

(11.5) 

From the master equation (II. 4) it follows that the matrix el- 
ements with respect to the basis |e), ji), jo), \g) of the state p 
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(11.6) 



Notice that if the initial state belongs to the class of density 
matrices of the form 
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(11.7) 



then p(f ) given by the solution of the equations (II. 6) also be- 
longs to that class. In this case one finds that 
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III. ENTANGLEMENT IN TWO-ATOMIC SYSTEMS 

A. Measure of entanglement 

In the case when subsystems of the total system are de- 
scribed by two-dimensional Hilbert spaces, the natural mea- 
sure of the amount of entanglement a given quantum state 
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contains i.e. the entanglement of formation | 



(III.l) 



where the minimum is taken over all possible decompositions 



Y.^kPk (111.2) 

k 



and 



E{P) = -tr [(tr^P) log2 (tr^P)] (III.3) 



can be analytically computed as a function of another quantity 
C(p) called concwrrence, which also can be taken as a mea- 
sure of entanglement C(p) is defined as follows 

C(p) = max (0,2p,^ax(p) -trp) (III.4) 

where /Jmax(p) denotes the maximal eigenvalue of p and 



with 



P = (p'/2ppl/2)l/2 



p = {^2 (E) 02) p {02(E) 02) 



(111.5) 



(III.6) 



The value of the number C(p) varies from for separable 
states, to 1 for maximally entangled pure states. Consider now 
the class (II. 7) of density matrices. With respect to the canon- 
ical basis fi, f2, fj, f4, the matrices (II. 7) have also the same 
form i.e. 



P = 



/pii pi4\ 

P22 p23 

P32 P33 

\P41 P44/ 



One can check that for this class 

C(p) = max{0,Ci,C2} 

where 



(in.7) 



Cl =2(|pi4| - VP22P33), C2==2(|P23|-VP11P44) 

(III.8) 



or 



Cl =2 Ip^gl - ^ (p„„ + p,,)2 - (p„, + p,„)2 



(111.9) 



C2 = V (p.« - PaaY- - {pm - P.u,)^ - 2 y/p^t 



when we use the matrix elements with respect to the collective 
basis. In the special case when p 14 = 0, Ci cannot be positive, 
so 

C = max{0,C2} 

Similarly, when P23 = 0, then 

C = max{0,Ci} 
B. Evolution of entanglement 

Suppose that the initial state of the two-atomic system be- 
longs to the class (II.7). Since the evolution given by the mas- 
ter equation (II. 4) leaves this class invariant, to compute en- 
tanglement at time / we can use formulas (III.9) and (II. 8). So 
we have 

C(p(f))=max{0,Ci(f),C2(f)} 

where 



and 



Ci{t)^2e-^^'\Peg\- 



-4^'(p«« + p..v-^) + ^ 



C2(f ) =e J (P.V.V - Paa)^ - (p«.v - P.va)^ - 



-e-4r'(p„, + p,„)2 



' 1 + e-8r'(- 1 + 2p,, + 2pg,y + Ae-^^' 



Pec + Pgg-7;-{9ee-Pgg)^ 



(HI. 10) 



(III.l 1) 



As an example, consider the following initial state 



The initial entanglement is equal to 



/O 0^ 

P22 P23 

p32 P33 

Vo oy 



(III. 12) 



C(p) = 2 IP23I = V (p.. - Paaf - {Pas " Psaf 



C{p{t)) = max{0,^.-2r'C(p) - - {\-e-'^')} (III.13) 
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From (III. 13) we see that there is the time fd(p) after which 
C(p(f )) becomes equal to 0. We may call fd(p) the time of dis- 
entanglement of a given initial state p. For initial state (III. 12) 



bounded from below. Now we can define the time of disen- 
tanglement rd(p) of a state p of two-atomic system as follows 



fd(p)-^ln 



C(p) + Jl+C(p)2 



(III. 14) 



As we show in the next section, the appearance of finite time 
of disentanglement for all initially entangled states, is the 
characteristic feature of the dynamics governed by the mas- 
ter equation (II.4). 



IV. TIME OF DISENTANGLEMENT 



fd(p) =inf5p 



(IV.3) 



Since the set of separable states is compact, Sp is closed and 
fd(p) S 5p. Therefore fd(p) may be also defined as the smallest 
time for which p(f) is separable. From the above discussion 
it is clear that for every initial state of the system, there exist 
finite time of disentanglement (which may be equal to for 
separable initial states). Consider now some explicite exam- 
ples. 



The evolution of states given by the semi-group generated 
by Lp^ has the following important properties: 

(i) it is local i.e. if p is separable, then p(f) is also separable 
for all t > 0, 

(ii) every initial state p evolves to a maximally mixed state 

Since maximally mixed state is separable and there is a neigh- 
bourhood of this state which contains only separable states, 
for an arbitrary state p the set 

Sp = {te [0,°o) : p(f) is separable} (IV.l) 

is always non-empty. Moreover, if 

Ep = {te [0,°o) ; p(r) is entangled} (IV.2) 

then Sp and Ep are disjoint and SpUEp = [0,°°). Notice that 
for every fi e Ep and every tj G Sp we have ti < tj, so Sp is 



1. Pure initial states. 

Since all pure entangled states ^ with concurrence C{^) = c 
are locally equivalent to the state given by the vector (see e.g. 

liol) 



\fl + 7f^, 0, 0, ^l-v/T^^ (IV.4) 

and our dynamics is local, it is enough to consider (IV.4) as 
initial state. The corresponding density matrix 



evolves into 








2 





1 



c 





(IV.5) 



,-4rf 



-iVT^e-^^' 





2ce 



-2Tt 





AT, 

1-e 



1 






-^^'-iVT^e-'-^'j 



(IV.6) 



One can easily check that 



(IV.7) 



This finite time of disentanglement should be contrasted with 
infinite time needed for decoherence process (see f5, U for 
discussion of different models). For pure initial states, we 
can introduce a decoherence rate 'k{P) which tells us how fast 
given pure initial state becomes mixed during the evolution 
L2J. If as a measure of mixedness we take the hnear entropy 

5iin(p) = l-trp2 



then 



m 



1 dSM{P{t)) 



2 dt 

In the case of state (IV.5) one obtains 



r=0 



HPi 



2r 



(IV.8) 



(IV.9) 



We see also that off-diagonal elements of P^ (f ) vanish asymp- 
totically with that rate. 

Another characteristic time of our evolution is connected with 
so called quantum non-locality i.e. the possibility of violat- 
ing Bell inequalities in quantum states 1 11^ il2il . It is known 
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that all pure entangled states violate Bell inequalities ifllll . but 
this is not true for mixed states. In the case of two two-level 
systems, there is an effective criterion for violating Bell in- 
equalities in mixed states QUI. For any density matrix p, 
take the real 3x3 matrix 

Tp = (fmn), f„m = tr(pa„(8)Om)n,»^= 1,2,3 



Define also real symmetric matrix 



(IV. 10) 



with eigenvalues mi , ui, ui,. Then p violates some Bell in- 
equality if and only if I14ll 



ot(p) > 1 



(IV. 11) 



and is always smaller then the time of disentanglement t(i{P<^) 
(see FIG. 1). 

We see that the evolution of pure initial states has the follow- 
ing remarkable properties: 

(a) in the time interval [O,fioc) the states (IV.6) are entangled 
and violate Bell inequalities, 

(b) for? S [fioc,fd), ^'<i>(f) are still entangled but do not violate 
any Bell inequality, 

(c) for all f > fd the states (IV.6) are separable, although deco- 
herence process takes infinite time. 

2. Some mixed initial states, 
(i) Consider the class of Werner states 



where 



m(p) — max {uj + Uk) 



(IV. 12) 



As a measure of nonlocality we may introduce a function 

n(p) max{0, hj(p) — 1} 
Initial state (IV.5) violates Bell inequality since 

m{P^) = l+c^ (IV.13) 

On the other hand, 

m{P^{t)) = e-^^' +c^e-^^' (IV14) 



decreases to 0, so there exists the time fioc after which (IV. 14) 
is smaller then 1 and nonlocality of initial state is lost. We see 
that this locality time for the state (IV.5) equals 



where 



W'± = (l-P)^+Pl^±)(^±l 



(IV. 16) 



V2 



[|O)0|O)±|1)0|1)] 



are maximally entangled pure states. It is known that W± are 
entangled for p > 1/3 and C(W±) — ^ • During the time 



(I V. 1 5 ) evolution W± become 



W±(f) = 



(l+pe 



\±2pe- 



-ATt 



2Tt 



I — pe 





-4n 



I — pe 




±2pe-^^" 

m 



(IV. 17) 



l+pe-^^'J 



and 



C(W±(0) = max <j Q,p [ e'^^' + i^e''^^' 



One finds that 



td{W±) 



2r 



In 



' + Vp{^+p) 



p> 



(IV. 18) 



(IV. 19) 



On the other hand, not all Werner states which are entangled, 
violate Bell inequalities. Nonlocal properties have only those 



states W± with p > -y= fl4ll . This nonlocality is lost when 
v2 

f >rioc(W±) 



(IV.20) 



But even in the interval -A= < p <l, this time is much smaller 
v2 

then time of disentanglement (FIG. 2) 
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FIG. 1: /(J and fjoc in units of [l/F] as functions of concurrence C, for 
pure initial states 
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mally entangled mixed states fvf], which have maximal value 
of entanglement for a given degree of inpurity measured by 
linear entropy 
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(g{c) Q c/2\ 
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\c/2 g{c)J 
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(IV.21) 
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FIG. 3: f(j and ?ioe in units of [l/F] as functions of c for maximally 
entangled mixed initial states 

where 



FIG. 2: f(j and in units of [1 /F] as functions of p for Werner initial 
states 



[1/3, ce [0,2/3] 
lc/2, ce [2/3,1] 



(ii) Similar computation can be done for the class of maxi- Direct calculations show that 



?d(p 



MEMSj 




^+2c2 + ^^(36c2+10)2-36 
1 - 2c + 4c2 + 2\/2cVl-2c + 2c2 



ce [0,2/3] 
ce [2/3,1] 



r 



(IV.22) 



One can also check that Pmems violates Bell inequality when 

c > — i= and for such values of initial entanglement, the local- 
v2 

ity time equals 



1 2 1 

floc(pMEMs) = ^ 1" 2c , > 



(IV.23) 



As in the previous cases, this time is always smaller then the 
time of disentanglement (FIG. 3). 



a noisy environment, modeled by the master equation (II.4) 
leads to the disentanglement of initially entangled states in fi- 
nite time. This time of disentanglement is (at least in examples 
considered above) the increasing function of initial entangle- 
ment - more entangled is the initial state, the longer period of 
time is needed to disentangle it. On the other hand, if the ini- 
tial entangled state can violate Bell inequalities, the period of 
time in which it still have nonlocal properties, is much shorter 
then the duration of the process of disentanglement. 



We have shown that interaction of two-atomic system with 
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